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Abstract 

The statistical properties of a Hamiltonian Hq perturbed by a localized scat- 
terer are considered. We prove that when Hq describes a bounded chaotic 
motion, the universal part of the spectral statistics are not changed by the 
perturbation. This is done first within the random matrix model. Then it is 
shown by semiclassical techniques that the result is due to a cancellation be- 
tween diagonal diffractive and off-diagonal periodic-diffractive contributions. 
The compensation is a very general phenomenon encoding the semiclassical 
content of the optical theorem. 
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In quantum systems, the chaotic or disordered nature of the classical motion is reflected 
in the statistical properties of the high lying eigenvalues and eigenvectors. For example, the 
spectral statistics of ballistic cavities are universal for energy ranges that are small compared 
to the inverse time of flight through the system. These universal properties are well described 
by random matrix theory (RMT) []I],[2|]. 

Consider now a perturbation imposed to a chaotic system. On the classical side, the 
dynamics is structurally stable, in the sense that a generic smooth perturbation leaves the 
dynamics chaotic. We are interested in the quantum mechanical effects of a particular 
class of perturbations that are non-classical. If the unperturbed motion is described by the 
Hamiltonian H acting in an N- dimensional Hilbert space, we consider Hamiltonians of the 
form 

H = H + XN\v)(v\ , (1) 

where \v) is a fixed vector. N is included in the perturbation for future convenience. The 
eigenvalues {u>i} of H satisfy the equation 

with {efc} the eigenvalues of H Q and Vk = (fk\v) the amplitudes of \v) in the eigenbasis of 

Rank-one perturbations like in Eq.(|]-0) appear in several contexts. The most common 
one is when a local short-range impurity or point scatterer is added to the system ||. 
The physical consequences of such a perturbation were studied for Fermi gases in 
the context of RMT || and for ballistic motion of particles in regular |7] and chaotic || 
cavities. Another context is in the physics of many body problems. In the basis of Hq, the 
perturbation simply reads Vkvf and is said to be separable. In a mean field approach, it is 
the simplest model leading to collective modes via the residual interaction ||. Although our 
results are general, we adopt for simplicity the language of the localized point scatterer. 

A local perturbation is purely wave-mechanical. For a system with / degrees of freedom, 
it represents a modification of the dynamics in a volume oc (2irh)f in phase space, that 

2 



tends to zero in the semi classical limit. For example, the addition of a point scatterer in a 
ballistic cavity leaves invariant the classical motion while quantum mechanically it induces 
wave effects like diffraction. The modifications of the eigenvalues by the perturbation are 
described by Eq.(^). A statistical analysis of the perturbed spectrum was done for the energy 
levels of a regular integrable cavity with a point scatterer in 0,|IIJ. It was demonstrated 



that a short range repulsion between the eigenvalues, different from RMT, is induced by the 
perturbation, thus considerably modifying the initial Poisson distribution. More recently, M. 
Sieber || has studied, using semiclassical techniques, the modifications by a point scatterer 
of the spectral statistics of chaotic systems. He showed that diffractive orbits produce 
finite contributions which may induce deviations with respect to the random matrix model. 
Whether this deviation really exists for chaotic systems, or on the contrary if there are 
other (non-diagonal) semiclassical contributions that cancel the purely diffractive terms is 
the question we answer here. 

We prove by two different approaches, namely a purely statistical model and a semiclas- 
sical calculation, that a local perturbation produces no deviations with respect to RMT. At 
the first place, assuming that the unperturbed eigenvalues and eigenvectors components in 
Eq.(^) are distributed according to RMT, i.e. their joint probability density are given by 
the standard formulae [T|Q] 

P{{e k }) oc II h - e/ , (3) 

i>j 

and 

//3iV\ 1-/3/2 

P (M) = II ( ^rj exp(-/5iV|^| 2 /2) , (4) 

we show that the joint probability density for the perturbed eigenvalues is exactly the same 
as the distribution of the unperturbed ones 

P{{u k }) ocJllut-Ujf . (5) 

i>j 

Here, (3 = 1 (resp. 2) for systems with (resp. without) time-reversal symmetry. In the 
second place, and to complete the analysis, a semiclassical calculation of the spectral form 



factor is considered. It is expressed as a double sum over all the periodic and diffractive 
orbits of the system. The latter are orbits that are scattered by the perturbation. In || the 
diagonal contribution of the diffractive orbits was obtained (cf Eq.(13) below). We compute 
the off-diagonal contribution coming from the interference of periodic and diffractive orbits, 
and find that this contribution exactly cancels the diagonal diffractive term. We thus recover 
the statistics of RMT. The basic physical ingredient at the basis of this cancellation is the 
unitarity of quantum scattering processes, i.e. conservation of the flux scattered by the 
impurity. Although our semiclassical result is less general than Eq.@ - it is only valid 
for the short time behavior of a two-point function -, it applies to systems with arbitrary 
diffraction coefficient not expressible as a Hamiltonian of the form ([I]). 

In chaotic and disordered system the local universal fluctuations of the spectrum are 
essentially related to the properties of the Jacobian (||). We ignore here problems related to 
the confinement of the eigenvalues, that are of minor importance for our purpose. We start 
the proof of Eq.([5]) from the joint distribution function of both the old and new eigenvalues, 
obtained in Ref. 

\ t i {jJ j\ 

with p = (3/2X. We restrict for simplicity to A > (A < is treated in the same manner). 
Eq.(^|) imposes the restrictions 6{ < u)i < (trapping). The distribution for the perturbed 
eigenvalues, uj iy is then defined as 

dei / de 2 .../ d e7 vP({e,},{^}) 
cxe-'S^n^-^) W (PiP) > (6) 

i>j 

with 

J i e pei dei r u N e peN de N 



and F(e) = flj \e— Ujl 1 ^ 13 ^ 2 . Using the standard expression for the Vandermonde determinant 
in W and integrating term by term one obtains 
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W(/3, p) = det 



r(*-i) 



i,j=l,...,N 



(7) 



where J^- = dllj is the i-i/i derivative with respect to p of 



I = I {0) 



> e pe de 



(8) 



For j = 1, uij-i = — oo. 

It is straightforward to check that the 1-,'s satisfy, for any j, the following differential 



equation [11 



2p 



This differential equation allows to write 



I J = 0. 



(9) 



with some coefficients a; t . W(f3,p) as defined in Eq.fl7|) is the Wronskian of this equation. It 
then follows that 



d p W = a N _, W , 

where in our case a/v-i = Y^i^i — flN/2p. Integration of Eq . (|T0"D leads to 



(10) 



W{P,p) 



p 



W/2 



ew(pJ2 



uO; 



When this result is replaced in Eq.(^) one gets 

pcm) <x Wo no*- 

i>j 
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3/ • 



(11) 



Wq is an integration factor, independent of p. We compute it from the asymptotic behavior 
of Ij when p — > oo. In this limit the integral in Eq.([|) may be evaluated explicitly 



lim L oc 



p— »oo 



p0/2 n^-k-^-r 



-/3/2 ' 



To leading order, 7^- = u;*/,,- and from Eq. (0) one gets 



\Ui — to. 



I/ 3 



WoocJIv 



From this equation and Eq. (|ll|) we recover the random matrix distribution function Eq. ([|) . 
A chaotic system coupled to the environment through a one-channel antenna has been 
considered in [|TIJ. The model is equivalent to Eq.(Q) but with imaginary A. For A —>■ oo 
the imaginary part of the perturbed energies is small and Eq.([5|) is obtained. Our method, 
which takes explicit care of the trapping problem, allows to prove this result for arbitrary A. 

In real physical systems agreement with random matrix theory is observed in a limited 
range. This universal behavior concerns correlations over energy ranges that are small 
compared to h/T min , with T min the typical period of the shortest periodic orbit. The above 
random matrix calculation establishes that the universal part of the spectrum is not changed 
by the presence of the scatterer. On the other hand, the non-universal behavior of the 
correlation functions occurring at scales of the order of, or larger than, h/T m[n are modified 
by the scattering center, since new diffractive orbits are introduced [p!3[|14[ . 

Let us now turn to a semiclassical treatment of the spectral correlations. These are based 
on trace formulae expansions of the density of states d{ui) = J2k ^(^ ~~ u k), written as a sum 
of smoothed plus oscillatory terms d = d + d^ osc ^ . We characterize the correlations by the 
spectral form factor defined as 

K(t ) = J°° ^ ^d (osc) (e + |) d^ - I)) exp fairjrd) . (12) 

The average indicated by brackets is taken over an energy window containing many quantum 
levels but whose size is small compared to E. We again consider a fully chaotic system with a 
point-like scatterer. In the geometrical theory of diffraction d^ osc ^ = <i^ osc ^ + d d ° sc \ where d p osc ^ 
and d d ° sc ^ are expressed as interferent sums over periodic and diffractive orbits, respectively 

MM- 

d ( ;f(E) =^^ex P (z^^-z|M ; (13) 

p,d 

with 

T 



27rft|det(M p - 1)| 1/2 
T d V{n,n') e-^+^ldetiVl 1 / 2 
4vr^(27r^)(/- 1 )/ 2 
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S Pj d(E) is the action of the periodic (resp. diffractive) orbits, T Pt d denotes their period, M p 
is the monodromy matrix of the periodic orbit, N is the matrix = d 2 Sd/dyi dyj (where y 
are coordinates orthogonal to the diffractive trajectory). P(n, ft') is the scattering amplitude 
of the scattering center located at xq with incoming n and outgoing ft' directions, defined in 
terms of the perturbed (G) and unperturbed (G ) Green's functions by the relation 

h 2 

G(x,x') = G {x,x') + —G (x,x )V(n, n')G (x , x') . 

Using the properties of the periodic orbits of chaotic systems, the diagonal contribution 
of S° sc ^ in Eq-dH) gives the short-time random matrix result K p (r) = (2//3)r [|15|]. The 



one-scattering contribution of the diffractive orbits in the same approximation is || 

with k the modulus of the wave vector at the impurity and a its total cross section 

a = J \V{n,n')\ 2 dQ dQ' . (15) 

(dJl is the solid angle element). For simplicity, we restrict the calculations to one scattering 
event (multiple scattering may be considered likewise). 

Our purpose is to compute the off-diagonal cross-term coming from the product of d p osc ^ 
and d d ° sc ^ in Eq.(^). The semiclassical expression for this contribution is 

K pd (r) = eMKS P - S d )/%) 5 (t - ?£±*k) + c.c.) . (16) 

p,d 

After energy smoothing, K p d has significant contributions only from orbits with close actions 
S p ~ Sd (having therefore approximately the same period). Pairs of orbits satisfying this 
condition may be constructed by considering the neighborhood of the forward scattering 
orbits. To each periodic orbit passing nearby the scatterer O we associate an "almost 
periodic" diffractive orbit that is similar to the periodic orbit but comes back to O with a 
slightly different momentum. In Eq.([TJ]) the double sum now involves all the possible pairs 
of trajectories constructed this way. Consider a surface of section that includes O and is 



transversal to the momentum of the periodic orbit when it comes nearby O. Let coordinates 
measured from O and momenta in the plane be denoted by (q,p). Consider all the periodic 
orbits of period T that cut the section through a differential element d^~ 1 q d?~ l p located at a 
distance q from O. The difference of action between these periodic orbits and the diffractive 
orbits associated to them as mentioned above is 

S P - S d = -(1/2) Qij q { q j , (17) 

with 

Qij = d 2 S/d gi dqj + d^/dq', d qj + d 2 S/dq[ dq'j, 

and q (q ') are initial (resp. final) coordinates on the surface of section. Moreover, one can 
show that 

| det Q\ = | det(M p - 1) det N\ cos 2 9 , (18) 

where 9 is the angle between the normal to the surface of section and the momentum of the 
diffractive orbit. 

By generalizing arguments used in the derivation of the Hannay - Ozorio de Almeida 
sum rule |l(| one can prove the following sum rule 

5(T-T P ) x(q P ,P P ) = /d^V^pxt^l nq] 
^ |det(M p -l)| S 1 ' 

where x(QjP) is a test function defined on the surface of section and (q p ,p p ) are the 

coordinates of the points at which the periodic orbit p crosses the surface of section. 

E = / d^xd^p5(i? — if(x, p)) is the total phase-space volume at energy E. From Eq. ([TED , 

using Eqs.(|17D and (0), we have 

J r 2 e i7r(/+l)/4 . 

'^ ^MH/J Vldet^-DHdetJVI 
xV*{n, fVje-^^^d^qd^p + c.c. 

Integrating the quadratic form in the exponent, taking into account Eq. (|IB|) , using the 
semiclassical density of states d = E/(27r^.)^ and the fact that the differential element for 
the momenta may be written d?~ 1 p = (hk)^' 1 cos9 dO, one obtains the final expression 
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2 / ; \ /-2 



fi[V*{n,n) -V(n,n))dfl . (20) 



2nf3 \2tt) 

This is the result for the cross-term contribution. Note that it depends only on T>(n,n), 
which translates the fact that in general interference terms between periodic and diffractive 
orbits can be large only for the diffraction in the forward direction. 

To make contact with Eq.flli]) we use a general relation valid for the elastic scattering 
on a finite range potential. The conservation of the flux scattered by the scattering center 
imposes a relation between the imaginary part of the scattering amplitude and the scattering 



cross section. This is the well known optical theorem [pTT| , that in / dimensions takes the 
form 

i[V*(n,n)-V(n,n)] = -^(J^j J |P(n,n')|W. 
Combining this relation with Eq. (^H) one gets our final result 

K pd (r) = -K d (r) . (21) 

The interference between periodic and diffractive orbits exactly cancels the diagonal contri- 
bution of the diffractive orbits, Eq. (JH|) . We recover from semiclassical methods, at least for 
a two-point function and short times, the RMT result. 

The two basic elements producing the cancellation are the sum rule (plj) and the optical 
theorem. Only the former is characteristic of chaotic systems, the latter being very general. 
The present semiclassical results may be extended by similar methods to multiple scattering 
events. In a wider context, it should be mentioned that this is one of the rare cases in which 
a calculation of off-diagonal contributions (whose role is essential in producing the correct 
result) is done explicitly for chaotic systems. 

We have concentrated on the fluctuation properties of eigenvalues of chaotic systems, 
and have demonstrated that they are unchanged by a local perturbation. This applies to 
high lying states, were the statistical hypotheses hold. On the opposite extreme, a local 
perturbation may lead to important modifications of the properties of the ground state of 
the system. Take for example a negative A. According to Eq.fl^), each perturbed eigenvalue 
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remains trapped by two unperturbed ones, except the ground state. The energy of the 
ground state may diminish arbitrarily with increasing |A| and, as can easily be shown, 
the associated wavefunction becomes more and more localized at the impurity. In our 
considerations we have ignored the presence of this "collective" mode. 

The authors are grateful for many useful discussions with O. Bohigas, M. Saraceno, M. 
Sieber and U. Smilansky. After completion of this manuscript we became aware of related 
semiclassical results obtained by M. Sieber. 
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